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Abstract 

Using integral equation theory of liquids to a binary mixed fluid lipid membrane, we study 
the membrane-mediated interactions between the macroions and the redistribution of neutral and 
charged lipids due to binding macroions. We find that when the concentration of binding macroions 
is infinitely dilute, the main contribution to the attractive potential between macroions is the line 
tension between neutral and charged lipids of the membrane, and the bridging effect also contributes 
to the attraction. As the relative concentration of charged lipids is increased, we observe a repulsive 
- attractive - repulsive potential transition due to the competition between the line tension of 
lipids and screened electrostatic macroion-macroion interactions. For the finite concentration of 
macroions, the main feature of the attraction is similar to the infinite dilution case. However, due 
to the interplay of formation of charged lipid - macroion complexes, the line tension of redistributed 
binary lipids induced by single macroion is lowered in this case, and the maximum of attractive 
potential will shift toward the higher values of the charged lipid concentration. 

Keywords: Two-dimensional fiuid membrane, membrane-mediated interaction, integral equation 
theory, line tension, bridging effect. 
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I. INTRODUCTION 



Generally, it is believed that, in the physiological condition, biomembranes are quasi two- 
dimensional (2D) fluid mixtures, composed of a wide variety of protein and lipid molecular 
species. The protein-protein, protein-lipid, and lipid-lipid interactions in and on the mem- 
branes are very important to the functions of the cell, such as ligand-receptor interactions 
(1), membrane rafts (2), and the formation of lipids domains and membrane budding (3), 
etc. Most of these problems are complex and not well understood yet. Recently, there is an 
increasing interest in understanding the electrostatic binding of charged macroions to the 
oppositely charged lipid membrane both experimentally and theoretically (4-13). Because 
of lipid's fluidity in the membrane, when the charged macroions are bound to the mem- 
brane, the oppositely charged lipids would migrate to the binding sites. This process would 
compete with the mixing entropic effects of different lipid species, and at equilibrium, the 
minimum of free energy is required. There are several theoretical works characterizing such 
a process (9-13). One of the models takes into account the entropy contributions of lipids 
and membrane-associated proteins to the free energy in the incompressible limit. Combined 
with the nonlinear Poisson-Boltzmann (PB) equation. May and Ben-Shaul have obtained 
the local lipid composition and the adsorption free energy in the single protein case, and ac- 
cessed the influence on these physical quantities when the concentration of binding proteins 
is flnite (9). By introducing a parameter x characterizing the extent of non-ideal lipid mixing 
in the mean field approximation, they were able to give the phase behavior of charged lipid 
membranes with the binding peripheral proteins. The existence of a two-dimensional phase 
separation suggests that there is an effective lateral attractive potential between the bind- 
ing proteins. However, the nonlinear PB theory itself does not introduce such a potential 
(14), and thus the authors argued that the attraction must be mediated by the membranes 
(10,12). To numerically solve the nonhnear PB equation simply and efficiently, most of the 
theoretical work mentioned above adopted cylindrical symmetry with the axis through the 
center of the macroion and normal to the charged membrane. Such a symmetric considera- 
tion make it difficult to calculate the effective pair potential between the macroions. In fact, 
there have been no systematic theoretical studies into detailed calculation of effective in- 
teractions between binding macroions. The specific form of membrane-mediated like-charge 
attractive potential has not been clarified, and the mechanism is still unclear. 

In the present paper, we establish a simple model of the fluid membrane and the binding 
macroions, to examine the lipid-mediated effective pair potential between binding macroions. 
Actually, the calculation of the effective pair potential is very important to understanding 
the phase behavior of soft matters. Recently, the effective potential in ionic solution that is 
beyond the description of the Derjaguin-Landau-Verwey-Overbeek (DLVO) theory receives 
great recognition (15), and it is believed that it will have further applications in biological 
systems (16,17). The most counterintuitive phenomenon is the effective attraction between 
like-charge objects. Most of these phenomena can be understood as the strong correlation 
effects of multivalent counter-ions (18-21). However, the long-range attraction between like- 
charge spheres conflned in parallel glass walls is still an open question both for experimental 
and theoretical studies (22-24). It is believed that such a long-range attraction is mediated 
by the substrate parallel glass walls because of confinement, but the mechanism is still not 
revealed theoretically (16). On the contrary, the theoretical understanding of effective poten- 
tial between the macroions mediated by the fluid membrane will be much more complicated, 
but is relatively unexplored. There are several approaches to get the effective potential theo- 
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retically. Here we will adopt the integral equation theory which is an accurate and powerful 
tool to explore the microscopic short-range structure of liquids (25,26), and recently, it has 
been applied successfully to explain the haloing effects of colloidal stabilization (27,28). 

II. MODEL AND THEORY 

It is well-known that biomembranes are a kind of lyotropic liquid crystals, and the lipid 

is an amphiphilic molecule composed of a hydrophilic head and one or two hydrophobic tail 
chains. In our model, wc take the membrane as two-dimensional fluids composed of neutral 
and charged lipid species. The interactions between the lipids can be separated into two 
parts which are head-head and tail-tail interactions. The head-head interactions are treated 
as 2D hard core, while the tail-tail interactions arc treated as 2D soft core to account for 
the repulsions between the tails originated from both enthalpy and entropy of chains. Thus 
the total pair potential Vai3{r) between a and (3 lipids with a lateral separation r can be 
written as: 

{oo, r<ro,h + rf3h, 

(^ap + Vap{r), r^h + ri3h<r < r^s + rps: (1) 

where the indices a and (= c or n), indicate charged or neutral lipids. Vah is a lipid's 
hard-core radius, and Vas is a lipid's soft-core radius, ea/3 is the soft repulsion between 
a and /3 lipids, therefore the interaction parameter characterizing the demixing extent of 
non-ideal mixed lipids is given by cUcn — ^cn ~ (^cc + ^nn)/'^- For the sake of simplicity, 
we set = e„„ = 0, and thus = ^cn- However, we should point out that the real 
interaction between chains may be more complex as presented in literatures (29-31). Va/3{r) 
is the effective potential between the lipids introduced by charged heads in the absence of 
binding macroions, and it is obvious that i'nn(r) = Vcn{r) — 0. 

The binding macroions can move, but confined in a two-dimensional plane which is par- 
allel to the membrane but with a separation h, as reported in literature (13). For simplicity, 
the macroion is modelled as a hard sphere, which is uniformly charged on the surface. The 
pair potential Vpp{r) between the macroions can be written as: 

^^^''> - \ Vpp{r), r > 2r„ 

where Tp is the hard-core radius of the macroion, Vpp{r) is the screened electrostatic inter- 
action in the absence of charged lipids. In the binding state, there is another pair potential 
to account for the interactions between the negatively charged lipids and positively charged 
macroions. Formally, it can be written as: 

Vpc{r) = Vpc{r,h), (3) 

where r is the lateral separation between the lipids and macroions, and h is the minimal 

membrane-macroion distance (sec Fig. 1). When h is fixed, the pair potential is only con- 
cerned with lateral distance r, so such a system can be effectively treated as a 2D fluid 
mixtures with pair potentials given above. Now, we can apply the two-dimensional inte- 
gral equation theory to our model. The Ornstein-Zernike (OZ) equations for homogeneous 
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mixtures in two-dimensional case are given by 

hafsir) = Ca/3{r) + ^cr^ Cau{\r - r'|)/i^;j(r')dr', (4) 

where hap{r) is defined by hapir) = gapif^) — 1- Qai^ir) is the radial distribution function, 
and hap{r) and Cap{r) are the total and direct correlation function between two particles 
of species a and /3, respectively, ai, is the number density of species u. The summation is 
over all species, and the integral is performed over the two-dimensional space. Eq.(4) can be 
closed by the well-known Percus-Yevick (PY) or hypernetted-chain (HNC) approximations or 
other thermodynamic self-consistent approximations (25,26,32) which have different bridge 
functions Bap{r). Formally, it can be written as follows: 

gal3{r) = Q-P^c8(r)+hc,i3{r)-Cc,i3{r)+B„i3{r) _ ^g-j 

In the HNC approximation, Bap{r) = 0, while in the PY approximation, Bap{r) = ln(l -|- 
^ap{r)) - lanir) where 7a/3(r) = hapir) - Cap{r). 

III. DETAILS OF CALCULATION 

The first step to do our calculation is to write down the specific form of pair potentials. 
Here, however, we must fall back on some approximations. The most important thing 
that we care is whether the charges are screened near the membrane in the physiological 
condition. One might argue that since the charges cannot penetrate into the membrane, 
the electro-static interactions between charged macroions or lipids near the membrane will 
not be screened. Such a statement does not take into account the specific properties of 
the membrane and the solution. We know that the hydrophobic part of the cell membrane 
has a dielectric constant Sm ~ 2.1 and a thickness around 6 ~ 7 nm, while the solution, 
in the physiological condition, has a dielectric constant Sg ~ 80 and the screening length 
A ~ 1 nm (33). When a charge Q occurs in the solution near the membrane (within 3 ~ 4 

£ — £ 

nm), it will form an image charge in the membrane with Q' = — —Q ^ Q. Now most 

of screening charges are near the membrane within 1 nm, and these charges would form 
images in the hydrophobic part of the membrane. Therefore the charges are surrounded by 
screening charges even on the membrane, and the electrostatic interactions must be more or 
less screened. This justifies the use of screened potentials. The screened potentials can be 
written as follows: 

+ [rp + hy 

B 



r 



where the Bjerrum length Lb — , A is the screening length of the solution, and 

4:TT£o£s 

f3 — l/ksT. Zc and Zp are the effective charges of charged hpids and macroions, respectively. 
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Substituting Eq.([6]) into Eqs.(l)-(3) gives the pair potentials we adopt here. Combined 
with the hard-core repulsion of macroions, we can study the effects of the screening length, 
charges, and macroion's size, etc. 

The three-component OZ equations are calculated in the k-space. With some algebra, it 
can be written in a form more suitable for numerical calculations, which is formally given 

by 

Kp{k) = fap{a,c{k)), (7) 

where c{k)= {ccn{k) , Cnp{k) , Cpdk) , Ccc{k) , Cnn{k) , Cpp{k)) and a= ((Tc, a„, (Tp). In the limiting 
(Tp 0, the OZ equations in the k-space reduce to a two-component OZ equations of lipids: 

^ ^f^^ ^ Cccjk) + tTn(5e«(fc) - 5cc(fc)Cn„(A:)) ^^^^ 
r _ Cnnjk) + ac{cln{k) — Ccc{k)Cnn{k)) 

""^ ' ~ D{k) ' ^ ' 

KM - If. (80 

and coupling equations with macroions: 

pel j = -^^^ , I aj 

7 /,N Cp„(A;)(l (^cCccik)) -\- (TcCcnCpc(yk) rriU\ 
ripn{K) - -^^^ , lyOj 

where D{k) = (1 - acCccik)){l - o-„c„„(A;)) - crcO-„c^„(A;). 

The two-dimensional Fourier transformation of the isotropic system can be defined by 

fir) = ^ r kfik)Mkr)dk, (10a) 

^TT Jo 

POO 

f{k) = 2tt rf{r)Jo{kr)dr, (10b) 



^0 

where Jo{x) is the zeroth-order Bessel function of the first kind. In numerical calculations, 
Eq. (fTUi) is truncated and discretized in r-space and /c-space. The discrete Fourier transfor- 
mation adopt here is a slow method, but satisfies the orthogonality conditions (34,35) which 
are important in our calculation. 

The OZ equations ([8]) are closed by the well-known PY approximation, which is good for 
short-range potential and predicts with a reasonable accuracy both structural and thermo- 
dynamic properties of the investigated system (25). For the short-range repulsion between 
the lipids studied here, it is obvious that the PY approximation is more suitable. On the 
other hand, macroion's effect can effectively be treated as an external field, and thus Eq. ([9]) 
is solved by combining with HNC approximation (26), which is better for the long-range 
potential. Finally, the integral equations are solved by using the Picard iteration scheme, 
and the equality can be rewritten as follows: 

clTir) = p[e-^^'^^^'''^+^-^^^^+^-^^''-^ - 1 - 7,^(r)] + (1 - p)<'^(r), (11) 

where p is an adjustable parameter to achieve the best performance of iterations. 
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IV. RESULTS AND DISCUSSION 



In the infinite-diluted case, the couphng between macroions has been ignored, and thus 
the distribution function of macroions can effectively be defined by gpp{r) = e~^^^p 
Because of the coupling between the macroions and the lipids, the explicit distribution 
function is written as follows: 

g^pl^r) = Q-l3Vpp{r)+ac J hpc{r')cpc{\r-r'\)dr'+an J hpn{r')cpn{\r-r'\)dr' ^ ^-|^2) 

namely the effective pair potential between macroions is 

PV^J^ir) = pVppir) -acj V(r')cpc(|r-r'|)dr' - an J V(r')cpn(|r-r'|)dr'. (13) 

Here, the effects of the lipids on the macroions have been included in the effective pair 
potential between the macroions. 

In the calculations presented below, we will restrict ourselves to the case that different 
species of lipids have the same size, — fnh — fcs = f^ns = 4r/3. The macroion's radius 

= 4r and the minimal mcmbrane-macroion distance h = 0.3r. The Fourier transformation 
is truncated at r ~ 90r. The choice of the soft-core range seems a bit arbitrary, however, 
since such a range is fixed during our calculation, it will not obscure the results given below. 
To keep in hue with the data appeared in Ref. (12), we choose Vch — Tnh — r — AA. 

In Fig. 2 we show the effective potentials (iV^^ {r) and the membrane- mediated potential 
l3Vpp{r) = PVpJ^{r) — f3Vpp{r) for varying values of the demixing strength f3ujcn between 
neutral and charged lipids in the physiological condition (A = 10 A, T = 300K). The lipid 
concentration is fixed to be acT^ = 0.055 for charged lipids and cr„r^ = 0.105 for neutral 
lipids. We see that with increasing the repulsion between the neutral and the charged lipids, 
the attraction between macroions is greatly enhanced. The result supports the conclusion 
made by Hinderliter (7,8) from experiments and Monte-Carlo simulations. Obviously the 
effective potential is composed of two parts: one is from the contribution of the screened 
potential between macroions, and the other is mediated by the substrate membrane. The 
membrane-mediated attraction between macroions is within 4 ~ 5 lipid's diameter from 
their surfaces. 

Figure 3 shows the total correlation functions hpc and hpn with the lateral separation 
between macroion and lipid for different values of Pucn- The correlations increases with the 
increase of PcUcn- Therefore, we can say that the attraction is mediated by the membrane 
through the correlations between the macroions and lipids. The resulting local density 
distribution (Ta{i^) of the a lipids which deviates from the averaged bulk density cTq,, can be 
obtained by the relation aa{r) = {hpa{r) -\- l)cra, when a macroion is in the center (r = 0). 
Interestingly, the density distribution of the lipids under the macroions is greatly modified 
by increasing the unmixing interaction fiujcn between neutral and charged lipids, while the 
electrostatic interactions between the macroions and charged lipids keep unchanged. As 
P'^cn = 0, the solid line in Fig. 3 shows that the enrichment of charged lipids appears below 
the macroion, just through the screened Coulomb potential. Such a process can be balanced 
by the reduction of the mixing entropy and the increase of the repulsions between charged 
lipids. As Pucn turns on, it has a tendency to drive the same lipids to aggregate together. 
This will greatly increase the macroion-lipid correlations, as shown in Fig. 3, by enhancing 
the aggregation of the same lipids, due to a pre-enrichment of charged lipids below the 
macroion at (5uJcn — 0. 
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The density profiles obtained here is similar to the results obtained from numerical solu- 
tions of the Poisson-Boltzmann equation (12,13). However in the present case, we have taken 
into account the correlation effects of species, and from the inset of Fig. 3(b), we find that 
after an enriched region of charged lipids, there is a depleted region of charged lipids where 
the number density of charged lipids is even smaller than the bulk case. Such a depleted 
effect of charged lipids will become stronger if the correlation effect is enhanced, as shown 
in Fig. 4 when Zc is increased to 2e. This is in agreement with the statement about the 
multivalent system where the charge correlations become important, and thus there exists 
some counter-intuitive effects (16), which are beyond the descriptions of the PB theory. 

For inhomogeneous mixing lipid system, if the interface between neutral and charged 
lipids is absolutely sharp on a molecular scale, the line tension will be directly proportional 
to the strength of unmixing (3uJcn- Therefore, to some extent. Figs. 2 and 3 indicate that such 
an attractive effect can be enhanced by the line tension between different lipids. However, the 
weak attraction exists, even for the case /ScUcn = where the line tension contributes little to 
the attraction. It indicates a deeper screening of the interactions between binding macroions 
in the presence of charged lipids. The origin of such a potential is due to the bridging effect 
of charged lipids, namely two macroions jointly attract to the same set of charged lipids 
(28). Intuitively, it seems that there exists another effect, the so-called depletion attractive 
interaction (36), to contribute to l3V^{r). If we take the charged lipids-rich clusters under 
the binding macroions as static, it is possible to introduce such a depletion effect by hard- 
core repulsions between the clusters and the surrounding lipids. However, such a picture is 
naive. The charged lipids-rich regime is not static, but quite dynamic. The charged and 
neutral lipids can move in and out of the regime quiet freely, although the density profiles 
remain unchanged. Therefore, there is no static interactions between the charged lipids-rich 
regimes and the surrounding lipids. Such a dynamic effect is similar to that of Karanikas 
and Louis (28). To illustrate it clearly, we eliminate the bridging effect by introducing a 
pseudo potential between the binding macroions and the charged lipids: 



Such a potential does not allow two macroions to attract the same set of charged lipids, but 
the charged lipids-rich regime still exists. When we apply such a potential to our model in 
the case (Sucn — 0, the weak attraction disappears. This means that the bridging effect is 
dominant and there is no depletion effect under the present dynamic picture. 

In Fig. 5, we change the relative composition of neutral and charged lipids, but the 
total number density of lipids is fixed. Interestingly, when is low, the effective potential 
is mainly repulsive, where the screened potential dominates over the interactions between 
macroions. With the increase of charged lipids, the membrane-mediated attraction becomes 
significant, and it competes with the screened potential at a distance about 2r - 6r and then 
forms an attractive well at that range. The attractive force saturates when acr"^ = 0.05. 
Further increase of charged hpids will weaken the attractive potential, because the screened 
potential begins to dominate the system again. 

To clarify the dominant contribution to the attractive potential, we can calculate the 
line tension by varying the number density of charged lipids, but jSuJcn is fixed. In the 
continuum limiting, the line tension is given by 7; = C{ucn) / dr(V?7c(*"))^, where ric{r) is 
the local composition of charged lipids, and is estimated by r)c{r) — adr) / (adr) + (7„(r)). 




(14) 
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C{u!cn) is a coefficient relevant to ujcn, here it is a constant, and thus we will present the 
result of 7/(r) in units of C{uJcn)- In Fig. 6, the solid square line tension curve 'yi{r)/C{uJcn) 
is plotted as a function of the number density of charged lipids for the infinite-diluted 
case (cTpT^ — > 0). The line tension first increases with cTcT^, and then is decreased. The 
maximum line tension almost correspond to the deepest attractive potential of lipid-mediated 
interactions between macroions in Fig.5, and they are of the same order. We can conclude 
that the attractive potential increases with the line tension, depending upon two factors: 
the demixing interaction and the relative component between neutral and charged lipids. 
On the other hand, the effective attractive potential between macroions is always enhanced 
by increasing charged lipid's concentration acT^, namely the attraction due to the bridge 
effect monotonously increases with cTcT^- Therefore, there is a cooperation between the line 
tension and the bridging effect which jointly assists the attraction between macroions at the 
beginning stage of the increase of UcT^, and for relatively large values of UcT^, a competition 
occurs, since line tension will decrease with further increase of charged lipids. The fact 
that the line tension behavior in Fig.6 is in accordance with the effective potential well 
vs (TcT^ in Fig.5, clearly indicates that the line tension plays a vital important role in the 
membrane-mediated potentials, as argued in the Ref. (12) in the continuum limits of the 
membrane. 

We now turn to discuss the effects of finite macroion concentration. In this case, the 
coupling between macroions cannot be ignored, and we have to solve the three-component 
OZ equations. The effective potential Vp^'^ between macroions can be obtained by an inverse 
process (37). After solving the integral equation Eq.(7), we can obtain the total correlation 
function hpp{r), and then the direct correlation function Ceff{k) can be solved from the 
reduced one-component OZ equation: 

ge//^^^ = T^^^TTT ■ (15) 
1 + aphpp{k) 

The effective potential can be solved from the HNC approximation (24,37): 

gppir) = e-^^-'('^)+'^-W-^^''('^). (16) 

As (Tp — s> 0, the present calculation reduces to the infinite-diluted case. 

In Fig. 6, we firstly show the line tension as a function of cTcT^ for different macroion 
concentrations. We notice that on increasing binding macroions, the line tension is gradually 
decreased and the position of the maximum line tension shifts toward the higher values of 
(Tc. These two behaviors are of the same origin, because there exists a competing interplay 
of different macroions to recruit charged lipids. The competition will partially suppress the 
modification of lipid composition profiles induced by a single macroion, and thus the line 
tension decreases, in comparison with the case apT^ — > 0. Furthermore, we can expected 
that, at lower macroion concentrations, the charged lipid's profile saturates under a single 
macroion and the line tension reaches the maximum at a certain value of acT^. However, if 
the macroion concentration is increased, the charged lipid's concentration will not saturate 
at a former value of iTcT^. To achieve the maximum line tension, we should further increase 
the charged lipid's concentration, i.e., the maximum line tension occurs with a higher value 
of (TcT^ with the increase of binding macroions. 

In Fig. 7 we depict the effective potential for varying values of acT^ when cxpr^ = 0.004. We 
can see that the change of effective potentials has the same tendency as in the infinite-diluted 
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case. When the concentration of charged hpids is low, the screened Coulomb potential 
dominates, and the effective potential is repulsive. In an intermediate concentration, the 
membrane-mediated attractive potentials compete with the screened Coulomb potentials 
and form an effective attractive well. With further increase of charged lipids, the mediated 
attractive potentials decay, and the screened Coulomb potentials dominate over the system 
again. If we calculate the line tension in this case, also, it has the same tendency as in the 
infinite-diluted case. However, we see from Fig. 7 that on increasing the binding macroions, 
the attractive potential saturates at a certain concentration with (TcT^ = 0.08, which deviates 
slightly from the value (JcT^ = 0.07 where the hne tension is maximum (Fig. 6). Such a 
deviation indicates that although line tension is still dominant for finite concentration of 
macroions, the bridge effect will become relatively stronger, compared to the infinite-diluted 
case. 



V. CONCLUSION 

In this paper, we have studied the membrane-mediated potential between the binding 
macroions which is induced by the interplay of the macroions and lipids. The process can 
be generally interpreted as follows. When the binding macroions are infinite dilute, the 
macroion will induce an enrichment of charged lipids under it. This process will be balanced 
by the decrease of the mixing entropy and the increase of repulsion between charged lipids. 
When there is a soft-core repulsion between charged and neutral lipids, the nonuniform 
distribution of lipids will introduce a line tension. The macroions prefer to aggregate together 
to reduce the line tension energy. This will compete with the screened Coulomb potential 
and form an attractive well in the effective potential. The bridging effect will also contribute 
to the attractive potential, but the line tension is dominant. Interestingly, as we increase 
the number density of charged lipids, but fix the total number density of lipids, we observe a 
repulsion - attraction - repulsion transition due to the competition between the line tension 
and screened electrostatic interactions. For the finite concentration case, we have to take 
into account the coupling effects between binding macroions. We find that with the increase 
of binding macroions, the line tension induced by a single binding macroion is gradually 
suppressed and the maximum of attractive potential will shift toward the higher values of 
charged lipid concentration Uc. The result can be attributed to the competition of recruiting 
charged lipids due to different binding macroions. The effective potential calculated in this 
case is not as faithful as the infinite dilution case. But in the case where the concentration 
of binding macroions is fixed, the result obtained is physically reasonable. 

This work was supported by the National Natural Science Foundation of China, No. 
10334020, No. 10021001, No. 20490220, and No. 10574061. 
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Figure captions: 

Figure 1: Schematic of a macroion bound to a binary fluid lipid membrane. The hpid is 
treated as a hard-sphere head with a radius of Tc/i (charged hpid) or rn/i (neutral lipid) and 
a soft-core repulsion tails. The binding macroions are treated as uniformly charged hard 
spheres with radius r^. There is a separation h from the membrane to the bottom of the 
macroion. The system is a three-dimensional system, but can be effectively mapped into a 
two-dimensional one when h is fixed. 

Figure 2: The effective pair potentials PVpJ-^{r) or PV^{r) (= PVp^'^{r)—j3Vpp{r)) change 
with Pojcn- The macroion's diameter is chosen to be 8r, and the effective charges of the 
macroions and charged lipids is Zp ~ 6e and Zc = — le, respectively. The mediated potential 
PV^{r) = j3VpJ^{r) — i3Vpp{r) becomes more attractive with the increase of PuJcn- This is a 
general behavior of the system. 

Figure 3: The total correlation functions hpnir) and hpdr) between the lipids and 
macroion for varying values of /9cJcn- In (b), the curves for the lateral distance r > 12r 
is enlarged, and rc-plotted in the inset, showing the correlation effects which will not appear 
in the mean field approximation. 

Figure 4: The total correlation function hpc{r) verse (3u!cn for Zc — — 2e. 

Figure 5: Effective pair potential (5V^^{r) verse the number density of charged lipids for 
fixed total lipids concentration (cTc + o"„)r^ = 0.16 and I3ujcn = 2.0. The attractive potential 
saturates at an intermediate concentration ((JcT^ = 0.05) of charged lipids. 

Figure 6: Single macroion induced line tension as a function of acr"^ for fixed total lipids 
concentration (uc -|- cr„)T^ = 0.16 and (5ujcn = 2.0. The solid square curve shows the case 
(TpT^ 0, and the line tension saturates at Uct"^ = 0.05. As the binding macroions are 
increased. The line tension decreases, and the position of maximum line tension moves 
toward the higher values of (JcT^. 

Figure 7: Effective pair potential PVpp{r) for varying concentrations of charged lipids 
with ((7c -I- (T„)t^ = 0.16, PcOcn — 2.0, and (7pT^ = 0.004. The attractive potential saturates 
at an intermediate concentration cTcT^ = 0.08. 
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